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25 Abstract 

26 Ubiquitin-like domains (Ubls) are now recognized as common elements adjacent to viral 

27 and cellular proteases; however, their function is unclear. Structural studies of the papain-like 

28 protease (PLP) domains of coronaviruses (CoVs) revealed an adjacent Ubl domain in Severe 

29 Acute Respiratory Syndrome CoV, Middle East Respiratory Syndrome CoV, and the murine 

30 CoV, mouse hepatitis virus (MHV). Here we tested the effect of altering the Ubl adjacent to 

31 PLP2 of MHV on enzyme activity, viral replication and pathogenesis. Using deletion and 

32 substitution approaches, we identified sites within the Ubl domain, residues 785-787 of 

33 nonstructural protein 3, which negatively affect protease activity, and valine residues 785 and 

34 787, which negatively affect deubiquitinating activity. Using reverse genetics, we engineered 

35 Ubl-mutant viruses and found that AM2 (V787S) and AM3 (V785S) viruses replicate efficiently 

36 at 37°C, but generate smaller plaques than WT virus, and AM2 is defective for replication at 

37 higher temperatures. To evaluate the effect of the mutation on protease activity, we purified WT 

38 and Ubl-mutant PLP2 and found that the proteases exhibit similar specific activities at 25°C. 

39 However, the thermal stability of the Ubl-mutant PLP2 was significantly reduced at 30°C 

40 thereby reducing the total enzymatic activity. To determine if the destabilizing mutation affects 

41 viral pathogenesis, we infected C57BL/6 mice with WT or AM2 and found that the mutant virus 

42 is highly attenuated, yet replicates sufficiently to elicit protective immunity. These studies 

43 revealed that modulating the Ubl domain adjacent to the PLP reduces protease stability and viral 

44 pathogenesis, revealing a novel approach to coronavirus attenuation. 

45 
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47 Importance 

48 Introducing mutations into a protein or virus can have either direct or indirect effects on 

49 function. We asked if changes in the Ubl domain, a conserved domain adjacent to the 

50 coronavirus papain-like protease, altered the viral protease activity or affected viral replication or 

51 pathogenesis. Our studies using purified wild-type and Ubl-mutant proteases revealed that 

52 mutations in the viral Ubl domain destabilize and inactivate the adjacent viral protease. 

53 Furthermore, we show that a CoV encoding the mutant Ubl domain is unable to replicate at high 

54 temperature or cause lethal disease in mice. Our results identify the coronavirus Ubl domain as a 

55 novel modulator of viral protease stability and reveal manipulating the Ubl domain as a new 

56 approach for attenuating coronavirus replication and pathogenesis. 

57 

58 Introduction 

59 Coronaviruses are emerging human pathogens. Severe Acute Respiratory Syndrome 

60 Coronavirus (SARS-CoV) caused the epidemic of 2002-2003, with ~10% case-fatality ratio (1). 

61 Middle East Respiratory Syndrome Coronavirus (MERS-CoV) is a pathogenic virus that was 

62 first identified in humans in 2012 (2). As of February 3, 2015 there have been 965 confirmed 

63 cases and 357 deaths (3). For SARS-CoV, the virus emerged from a reservoir in bats, replicated 

64 in an intermediate host (civet cats), and spread to humans. The epidemic strain of SARS-CoV 

65 evolved for efficient human-to-human spread (4-6). Public health measures of isolation of 

66 infected individuals led to the cessation of the epidemic in humans; however, SARS-like viruses 

67 remain in bat reservoirs (7-9). For MERS-CoV, dromedary camels are now suspected as the 

68 likely zoonotic source for transmission to humans since MERS-CoV sequences with 99% 

69 nucleotide identity to human MERS-CoV isolates have been detected in respiratory samples 
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70 from camels (10). Although there are reports of human-to-human transmission of MERS-CoV 

71 (11,12), current strains seem to cause mostly lower respiratory tract disease and are not as highly 

72 transmissible as SARS-CoV (13). Other human coronaviruses (HCoV-229E, HCoV-OC43, 

73 HCoV-NL63, HCoV-HKUl) are endemic in the human population and are the causative agents 

74 of upper and lower respiratory tract disease and croup (14-17). To date, there are no FDA 

75 approved antiviral drugs or vaccines to fight human coronavirus-induced disease. In addition, 

76 the potential exists for coronaviruses to emerge into the human population from endemic 

77 reservoirs in bats or other animals (18). Identifying viral components critical for efficient 

78 replication and manifestation of disease will facilitate antiviral drug and vaccine development. 

79 Coronaviruses are RNA viruses that encode a replicase polyprotein at the 5’-end of the 

80 positive-strand genome. Upon virus entry, the genomic RNA is translated to produce replicase 

81 polyproteins (ppla and pplab), which are processed by virally-encoded proteases. Depending on 

82 the virus species, the ppla encodes one or two papain-like proteases (PLPs) and one 3C-like 

83 proteinase (3CLpro or Mpro). These proteases process the replicase polyproteins into non- 

84 structural proteins (nsps) that assemble with cellular membranes and facilitate virus replication 

85 (reviewed in (19)). Mouse hepatitis virus (MHV), the murine coronavirus used in this report, is 

86 commonly used as a model system to study the replication and pathogenesis of coronaviruses 

87 (20). 

88 The MHV replicase product nsp3 consists of multiple domains including two PLP 

89 domains (PLP1 and PLP2), two predicted ubiquitin-like domains (Ubl), an acidic region (Ac), 

90 ADP-ribose-l”-phosphatase (ADRP), nucleic acid-binding domain (NAB), coronavirus group 2 

91 marker domain (G2M), transmembrane segment (TM) and coronavirus highly conserved domain 

92 (Y) (reviewed in (21), shown in Fig. 1A). The Ubl-1 domain has been shown to interact with 
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93 nucleocapsid (N) protein, which is important for virus replication (21, 22). Previous studies 

94 revealed that MHV PLP1 activity is required for processing the polyprotein at the nspl/nsp2 and 

95 nsp2/nsp3 sites (23). Further, the catalytic activity of PLP1 is required for efficient virus 

96 replication (24). PLP2 was shown to recognize and process a LXGG motif and cleave the 

97 replicase polyprotein at the nsp3/nsp4 junction (25). The LXGG recognition site is similar to the 

98 RLRGG recognition site of cellular deubiquitinating enzymes. Lindner et al., 2005 was the first 

99 to predict that CoV PLPs could be multifunctional enzymes with both protease and 

100 deubiquitinating (DUB) activity (26). Indeed, further studies revealed that SARS-CoV PLpro, 

101 MERS-CoV PLpro, HCoV-NL63 PLP2, PEDV PLP2 and MHV PLP2 are multifunctional 

102 enzymes exhibiting protease, DUB and delSGylating (ability to deconjugate interferon 

103 stimulated gene 15 protein, ISG15 from substrates (27-34)) activities. CoV PLP activity is 

104 required for processing the replicase polyprotein and predicted to modulate the innate immune 

105 response by deubiquitination of signaling molecules activated by pattern recognition receptors 

106 such as RIG-I and MDA5 (29, 31, 35, 36). Structural studies of the papain-like protease domain 

107 of SARS-CoV (37, 38) and more recently of MERS-CoV (39) and MHV (Chen et al., in 

108 preparation), show that these enzymes belong to the ubiquitin-specific proteins (USP) family of 

109 deubiquitinating enzymes and reveal the presence of a ubiquitin-like domain (Ubl) located 

110 upstream of the protease domain (38). However, the role of the Ubl domain in modulating the 

111 enzyme activity of MHV PLP2 is unknown. 

112 Here, we investigated the role of the Ubl domain (designated Ubl-2) adjacent to PLP2 for 

113 its function in MHV replication and pathogenesis. We generated proteases and viruses 

114 containing mutations within the Ubl domain and found that these mutations decreased PLP2 

115 activity and stability. Further, we found that these mutations resulted in decreased virus 
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116 replication and a marked attenuation of virulence. Immunization with the Ubl mutant virus 

117 protected mice against challenge with wild-type virus. Overall, our data demonstrate for the first 

118 time that manipulation of the Ubl domain adjacent to a viral protease can ameliorate viral 

119 pathogenicity in vivo. 

120 

121 Materials and methods 

122 Cells. HEK293T cells were cultured in Dulbecco's Modified Eagle Medium (DMEM) with 10% 

123 fetal bovine serum (FBS) and 2% L-glutamine. DBT cells were cultured in Minimal Essential 

124 Medium with 5% FBS, 2% L-glutamine, and 10% Tryptose Phosphate Broth (TPB). BHK-R 

125 cells were kindly provided by Mark Denison (Vanderbilt University Medical Center) and 

126 maintained in DMEM media supplemented with 10% FBS, 2% L-glutamine, and 0.8mg/ml 

127 G418. 

128 Plasmids and mutagenesis. PLP2 (amino acids 1525-1911) sequence in frame with a V5 

129 epitope tag was codon-optimized and synthesized by Gene Script (Piscataway, NJ) (sequence 

130 available upon request). Codon-optimized synthetic PLP2 sequence was cloned into pCAGGS- 

131 MCS vector. For mutagenesis, an overlapping PCR strategy was used with primers described in 

132 Table 1. The introduced mutations were verified by sequencing. The nsp2/3-GFP expression 

133 plasmid was kindly provided Ralph Baric (University of North Carolina). The Flag-Ub plasmid 

134 was kindly provided by Adriano Marchese (Loyola University Chicago). 

135 Protease and Deubiquitinase (DUB) Activity Assays. To determine catalytic activity of the 

136 PLP2 constructs, 70% confluent HEK293T cells in 12-well CellBIND plates (Corning) were 

137 transfected using ThmsIT-LTl Reagent (Mims) according to the manufacturer’s protocol. For 
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138 the protease assay, the cells were transfected with 25ng nsp2/3-GFP plasmid and 300ng 

139 pCAGGS-PLP2-V5 expression plasmids (wild-type and various PLP2 mutants). To assess DUB 

140 activity the cells were transfected with 600ng Flag-Ub plasmid, and pCAGGS-PLP2-V5 

141 expression plasmids (wild-type and various PLP2 mutants). At 24 hours post-transfection cells 

142 were lysed with 300pL of lysis buffer A containing 4% SDS and 3% dithiothreitol (DTT). 

143 Proteins were separated by SDS-PAGE and transferred to PVDF membrane in transfer buffer 

144 (0.025M Tris, 0.192M glycine, 20% methanol) for 1 h at 55V at 4°C. Following this, the 

145 membrane was blocked using 5% dried skim milk in TBST buffer (0.9% NaCl, lOmM Tris-FICl, 

146 pFI7.5, 0.1% Tween 20) overnight at 4°C. The membrane was incubated with polyclonal rabbit 

147 anti-GFP antibody (Life Technologies) at a dilution of 1:2000 for the protease assay, or mouse 

148 anti-flag (Sigma) at the dilution of 1:2000 for the DUB assay. The membrane was washed 3 

149 times for 15 minutes in TBST buffer followed by incubation with secondary donkey-anti-rabbit- 

150 FIRP antibody at a dilution of 1:2000 (Amersham) for protease assay, or goat anti-mouse-FIRP 

151 antibody at a dilution of 1:5000 (Amersham). Then the membrane was washed 3 times for 15 

152 minutes in TBST buffer. Detection was performed using Western Lighting Chemiluminescence 

153 Reagent Plus (PerkinElmer) and visualized using a FluoroChemE Imager (Protein Simple). To 

154 verify expression of the PLP2 constructs, membranes were probed with mouse anti-V5 (Life 

155 Technologies) antibody at the dilution 1:5000. Mouse anti-calnexin (Cell Signal) antibody at a 

156 dilution 1:2000 was used as a loading standard. 

157 Biosensor Live Cell Assay. To determine protease activity of the Ubl mutants, the previously 

158 described protocol was used (40). Briefly, F1EK293T were transfected with 37.5ng pGlo- 

159 RLKGG construct and 50ng of PLP2 expression plasmids. At 18 hours post-transfection, 
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160 GloSensor (Promega) reagent diluted 1:50 in DMEM (10% FBS) was added. The luminescence 

161 was measured using a luminometer (Veritas) every hour for 5 hours. 

162 Generating Ubl-2 mutant viruses. To introduce mutations into MHV A59, we used a 

163 previously described method (40). Briefly, plasmid encoding the MHV B subclone was 

164 mutagenized using primers described in table 1. Plasmids encoding the complete virus genome 

165 were digested with restriction enzymes, gel purified, and ligated using T4 ligase at 16°C 

166 overnight. The ligation reaction was isopropanol precipitated and in vitro RNA transcription was 

167 performed using a mMESSAGE mMACHINE Kit (Ambion) according to the following 

168 protocol: 40.5°C for 25 min, 37.5°C for 50 min, 40.5°C for 25 min. RNA was electroporated into 

169 BHK-R cells, and the electroporated cells were seeded onto DBT cells. The supernatant was 

170 harvested 36 hours post-electroporation and plaque assay was performed as described previously 

171 (24). RNA was extracted from infected DBT cells at 8 hours post-infection using RNeasy Mini 

172 (Qiagen) and cDNA was generated using RT 2 First Strand Kit (Qiagen). PCR was performed 

173 using replicase primers (Table 1) and purified PCR product was sequenced (amino acids 747- 

174 848). AM2 and AM3 were plaque purified and AM2 was subjected to deep sequencing. 

175 Temperature Shift Experiment. DBT cells in 6-well plates were infected with 0.1 MOI of 

176 wild-type MHV or AM2 at 37°C. At 2, 4, or 6 hours post infection, cells were moved to 39.5°C 

177 until 14 hours post-infection when supernatant was harvested and virus titer was determined by 

178 plaque assay at 37°C as described previously (25). 

179 Expression and purification of wild-type and V787S mutant proteins. The PLP2 sequence 

180 was cloned from pCAGGS-MCS-PLP2 (described above) into LIC vector pEV-L8, which is a 

181 modified pET-30 plasmid. The expression of wild-type PLP2 and V787S mutant was performed 
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182 using Escherichia coli strain BL21 (DE3). Cultures were grown in LB medium supplemented 

183 with kanamycin (50 ug/ml) at 37°C until the optical density at 600nm (A600) reached 0.6. PLP2 

184 expression was then induced with O.lmM isopropyl-P-D-thiogalactopyranoside (IPTG) at 25°C 

185 for 6h. Wild-type PLP2 was induced at 25°C for 6 h, while the V787S mutant was induced at 

186 18°C overnight. Cells were harvested by centrifugation at 4,690 x g for 20 min at 4°C and stored 

187 at -80°C until use. Cell pellets from 1 L culture were then resuspended in 40mL buffer A (25mM 

188 Tris pH 7.0, 500mM NaCl, 20mM imidazole, 5mM PME) supplemented with dissolved flakes of 

189 lysozyme and DNase, lysed through sonication and centrifuged at 28,880 x g for 30 min (4°C). 

190 The supernatant was filtered through a 0.45pm membrane (Millipore) and loaded onto a 5ml Ni 

191 HiTrap HP column (GE healthcare) pre-equilibrated with buffer A. Then the column was washed 

192 with buffer A supplemented with 5% buffer B (25mM Tris pH 7.0, 500mM NaCl, 500mM 

193 imidazole, 5mM PME) until the UV was back to the baseline. The protein was eluted through a 

194 gradient of 5%-100% buffer B in 30 column volumes (CV). Fractions were collected and then 

195 pooled after enzymatic activity and purity assessment. The n-terminal (His)8-tag was then 

196 removed by TEV protease (His-tagged) cleavage by incubating PLP2 and TEV-protease together 

197 overnight at 4°C while dialyzing into buffer C (25mM Tris pH 7.0, lOOmM NaCl, lOmM PME). 

198 Then, free His-tag and TEV protease and uncleaved PLP2 were then separated from cleaved 

199 PLP2 by running the sample over a Ni 2 -charged HiTrap column. The flowthrough was 

200 collected and then concentrated using Millipore Micron concentrators to a volume of less than 2 

201 mL. The concentrated sample was loaded onto a Superdex-75 Hiload 26/60 column (GE 

202 Healthcare) pre-equilibrated with buffer D (50mM HEPES pH 7.0, lOOmM NaCl, lOmM DTT), 

203 and eluted at a flow rate of 2 ml/min. Fractions containing active enzyme at high purity, as 
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204 judged by SDS-PAGE, were pooled, concentrated and flash-frozen in 2% glycerol using liquid 

205 nitrogen for storage at -80°C. 

206 Temperature inactivation of wild-type and V787S purified protein. Wild-type PLP2 and the 

207 V787S mutant protein were incubated at 25°C for different time periods (0-50 min). At each time 

208 point, the specific activity of both enzymes was measured at 25°C using a Synergy Multi-Mode 

209 Microplate Reader (BioTek) with 50 uM RLRGG-AMC as the substrate and 3 pM of each 

210 enzyme. The assay buffer used contained 50mM HEPES, pH 7.0, 0.1 mg/ml bovine serum 

211 albumin (BAS) and 2mM DTT. The experiments were performed in triplicate in a final volume 

212 of 100 pi using the 96-well Corning Costar black microplates. Similar experiments were carried 

213 out when the enzymes were incubated at 30°C. To analyze the kinetic data, the ratio of the 

214 reaction rate at time=t to the rate at time=0 was plotted on a logarithmic scale against incubation 

215 time. Kinetic data of the V787S PLP2 incubated at 30°C were fitted to a first-order exponential 

216 decay model (Rate/Rateo = e" kt ), from which the inactivation rate constant k,„ ac( and half-life ti /2 

217 were determined. The wild-type data were fit to a line since no significant temperature 

218 inactivation was observed. 

219 Thermal melting temperature (T m ) analysis using circular dichroism (CD). Thermal melting 

220 analyses of the wild-type PLP2 and V787S mutant was carried out with a Chirascan circular 

221 dichrosim (CD) spectrometer (Applied Photophysics) equipped with a temperature control 

222 system (Quantum Northwest Inc.) by monitoring the CD signal at 220nm while increasing the 

223 temperature at a step interval of 0.4°C and at a rate 0.5°C/min. Two ml of Protein samples at 

224 lpM in buffer with 0.1M potassium phosphate (pH 7.5) was contained in a 10 mm quartz cell 

225 (Starna Cells) with magnetic stirring. Thermal scans were performed in three independent 
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226 experiments for both wild-type and V787S mutant MHV PLP2. The melting temperatures (T m ) 

227 were calculated as the first derivative peak using the program SigmaPlot. 

228 Generation of bone-marrow derived macrophages and virus infection. Femur and tibia from 

229 8-week old C57BL/6 female mice were prepared and bone-marrow was isolated, and cultured in 

230 IMDM with 10% FCS, 1% pen-strep, 0.1% beta-mercaptoethanol and 5% MCSF obtained from 

231 L929-supernatants. Additional medium was added at day 3 post-isolation and cells were seeded 

232 at day 6 at a density of 2 * 10 5 cells/well. The following day, cells were infected with wild-type 

233 icMFlV A59 and AM2 at an MOI of 1. Medium was changed after two hours and supernatant 

234 was harvested at various time points post-infection for titration. Supernatant was titrated on 

235 murine fibroblast L929 cells (in MEM containing 10% FCS and 1% pen-strep). 

236 Mouse studies. C57BL/6 mice were purchased from the Jackson Laboratory. The mice were 

237 maintained, and experiments performed at Loyola University Chicago in accordance with all 

238 federal and university guidelines. 4-week old C57BL/6 male mice were anesthetized with 

239 ketamine/xylazine prior to intracranial injection with 600 PFU of wild-type icMHV A59 or AM2 

240 mutant. Weight loss and survival were monitored on a daily basis. Mice were sacrificed when 

241 they lost 25% of their initial body weight. In other experiments, mice infected with AM2 at 4 

242 weeks of age, and age-matched controls were challenged with 6,000 PFU of wild-type icMHV- 

243 A59 at 9 weeks post primary infection. Survival and weight loss were monitored daily. Mice 

244 were sacrificed when they lost 25% of the initial body weight. For histology livers were fixed in 

245 4% formalin and paraffin embedded. Sections were stained with hematoxylin and eosin, and 

246 analyzed by microscopy. 
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247 qRT-PCR analysis. Brains and livers of 4-week old male mice infected with wild-type icMHV 

248 or AM2 were harvested and homogenized at day 3 and 5 post infection and RNA was extracted 

249 using RNeasy Mini (Qiagen). Reverse transcription was performed using 2 ug of total RNA and 

250 the RT 2 First Strand Kit (Qiagen) according to manufacturer’s protocol. 1 pi of cDNA was used 

251 to set up qRT-PCR reaction according to the manufacturer’s protocol using Single Primer Assay 

252 for IL-6 (SABiosciences), or replicase primers (table 1). Ct values were normalized to the 

253 housekeeping gene (Actin). 

254 

255 Results 

256 Identifying the ubiquitin-like domain 2 (Ubl-2) as a modulator of MHV papain-like 

257 protease activity 

258 To identify domains and residues within nsp3 that modulate protease activity, we 

259 generated a synthetic, codon-optimized parent construct encoding the protease (PLP2, residues 

260 699-1087) in-frame with a V5 epitope tag (Fig. 1A) and introduced a series of mutations into the 

261 predicted ubiquitin-like domain 2 (Fig. IB). Protease activity was evaluated using an established 

262 trans-cleavage assay (Fig. 1C). Our goal was to identify a domain that modulated, but did not 

263 inactivate protease activity. As a negative control, we generated a PLP2 catalytic mutant 

264 (C890A) that was inactive in the trans-cleavage assay. We generated a quadruple mutant (aa 

265 785-788 VDVL/SSSS) and specific single mutants: V785S, D786K, and V787S. We found that 

266 substitution of specific residues within the ubiquitin-like domain 2 adjacent to PLP2 consistently 

267 affected protease activity and reduced the abundance of cleavage product detected in the trans- 

268 cleavage assay when compared to the activity of PLP2-WT (Fig. 1C). To determine the degree to 

269 which mutants’ protease activity declined we used the recently described GloSensor assay (41). 
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270 We found that, consistent with trans-cleavage assay results, all Ubl-2 mutants had decreased 

271 protease activity. However, VDVL/SSSS, V785S, and V787S showed the highest defect in 

272 protease activity, while D786K showed intermediate phenotype (Fig. ID). In these experiments 

273 similar levels of PLP2 wild-type and mutants were expressed in each sample. These results 

274 implicated the Ubl-2 domain as a modulator of viral protease activity in MHV. Further, we 

275 determined the deubiquitinating (DUB) activity of the Ubl-2 mutants. We found that quadruple 

276 mutant and the single mutants V785S, and V787S had reduced DUB activity compared to wild- 

277 type PLP2. In contrast, D786K mutant had similar DUB activity to wild-type PLP2 (Fig. IE). 

278 Therefore, we focused on the Ubl-2 V785S and V787S mutants since these Ubl mutants 

279 exhibited the most significant defects in the cell-based assays. 

280 Ubl-2 mutant virus MHV-AM2 is temperature sensitive 

281 To determine if substitutions in the Ubl-2 domain affected virus replication, we used 

282 reverse genetics (42) to introduce specific mutations into the viral genome and assessed viral 

283 replication. Since the phenotype of the mutations in the PLP2 assay (Fig. 1) can be different 

284 from the phenotype of the mutant in the context of replicating virus we generated viruses 

285 encoding a quadruple mutation VDVL/SSSS, designated MHV-AM1, a single mutation V787S, 

286 designated MHV-AM2, and another single mutation V785S, designated MHV-AM3 (Fig. 2A). 

287 The presence of the mutations was confirmed by sequencing within the PLP2 domain. MHV- 

288 AMI was highly impaired, generated pinpoint size plaques, and grew to low titer (lxl O' 

289 PFU/ml). The low titer and pinpoint plaques generated by AMI are likely due to the quadruple 

290 mutation in PLP2; however, we cannot rule out the possibility that additional mutations outside 

291 of PLP2 contribute to this phenotype. 
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292 In contrast, AM2 and AM3 viruses with single amino acid substitutions within the Ubl-2 

293 domain replicated to titers similar to the parent infectious clone. We found that both AM2 and 

294 AM3 viruses generated slightly smaller plaques than those generated by wild-type MHV A59. 

295 The small plaque phenotype can be associated with temperature sensitivity, so we evaluated 

296 relative plaque size of wild-type and AM2 by performing plaque assays at 39.5°C and 40.5°C. 

297 We found that AM2 produced small plaques at 39.5°C (about 50% reduction compared to wild- 

298 type), but was unable to generate plaques at 40.5°C (Fig. 2B), consistent with a temperature 

299 sensitive defect. AM3 also produced slightly smaller than wild-type plaques but bigger than 

300 AM2 plaques at 39.5°C (about 20% reduction compared to wild-type) (data not shown). Because 

301 the temperature sensitive phenotype was more severe for AM2 and because V787 is conserved in 

302 other betacoronaviruses such as SARS-CoV and MERS-CoV (Fig. IB), we focused our studies 

303 on AM2. We verified the engineered sequence of AM2 by deep sequencing the viral genomic 

304 RNA. The 2-nucleotide substitution (gtc to ccc) of V787 was confirmed, and no additional 

305 changes were detected (GenBank # in progress). 

306 To determine the replication kinetics of AM2, we infected DBT cells with isogenic wild- 

307 type MHV-A59 or AM2 virus at MOI of 0.1. The cells were incubated at 37°C, supernatant was 

308 harvested and virus titer was determined by plaque assay. We found that at 37°C, AM2 

309 replicated with similar kinetics to wild-type virus (Fig. 3A). In addition, no replication defect 

310 was observed in bone marrow derived macrophages, important primary target cells during MHV 

311 infection (Fig. 3B). However, if infected DBT cells were shifted to a higher temperature (39.5°C) 

312 there was a progressive reduction in viral titer with increased time at the higher temperature (Fig. 

313 3C). When the cells were incubated at 37°C for 6 hours and then shifted to the higher 

314 temperature, the difference between the titers of WT and AM2 was about half a log. In contrast, 
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315 when the infected cells were incubated at 37°C for two hours and then shifted to 39.5°C, the titer 

316 of AM2 was two logs lower compared to wild-type virus. Taken together, these data show that 

317 AM2 encoding V787S has a temperature sensitive defect in viral replication. 

318 Nsp3-V787S substitution impairs PLP2 thermal stability 

319 To investigate the temperature sensitive phenotype conferred by the V787S substitution, 

320 we expressed the wild-type and V787S mutant of MHV PLP2 in E. coli and then purified the 

321 enzymes to homogeneity using a multi-step chromatography procedure (Fig. 4A). The thermal 

322 stability of the wild-type and V787S purified enzymes was then determined at 25°C and at 30°C 

323 by measuring the specific activity of the enzymes as a function of incubation time using the 

324 peptide-RLRGG-AMC as a substrate. We found that when wild-type and V787S mutant PLP2s 

325 were incubated at 25°C, both enzymes were stable over the incubation time course and had 

326 similar specific activity (Fig. 4B). The activities in Fig. 4B have been normalized to 100%. In 

327 contrast, when wild-type and V787S mutant PLP2s were incubated at a higher temperature 

328 (30°C), the specific activity of V787S PLP2 decreased rapidly over time whereas the wild-type 

329 protein maintained full activity (Fig. 4B). These results indicate that the V787S enzyme is 

330 thermally unstable at 30°C with a half-life of27.7±l.l min. 

331 To determine the relative structural stability of the purified proteins over a wider 

332 temperature range, we performed thermal-denaturation experiments using circular dichroism 

333 (CD) analysis to determine the melting temperatures (T m ) of each enzyme. Consistent with 

334 Figure 4B, the V787S mutant exhibited diminished thermostability (Fig. 4C). The determined 

335 melting temperature was 6.8°C lower for the V787S mutant protein than for the wild-type PLP2, 

336 suggesting that the V787S mutant PLP2 is less stable and unfolds at a lower temperature than 

337 wild-type PLP2. Taken together, our data indicate that the V787S substitution diminished the 
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338 overall protein thermal stability resulting in decreased enzymatic activity at elevated 

339 temperatures. 

340 AM2 exhibits reduced pathogenesis and elicits a protective response in mice 

341 To determine the effect of the Ubl-2 domain mutation on virus pathogenesis, we infected 

342 4-week old male C57B/L6 mice with 600 PFU of wild-type or AM2 virus by intracranial 

343 injection. We monitored weight loss over the course of infection and mice were humanely 

344 euthanized if they lost greater than 25% of initial body weight. We observed that while 100% of 

345 mice infected with wild-type virus succumbed to infection by day 7, all mice infected with AM2 

346 survived the infection (Fig. 5A). We confirmed that the mutation in AM2 was maintained upon 

347 mice infection by sequencing within the PLP2 domain of recovered AM2. 

348 We detected similar viral titers in the brains of wild-type or AM2-infected mice at day 3 

349 or day 5 post infection (Fig. 5B). In contrast we found levels of viral replicase and IL-6 were 

350 reduced in the livers of AM2 infected mice at day 5 post infection compared to wild-type MF1V 

351 infected mice (Fig. 5C). There was no difference in CCL2 or IFN-a (data not shown). 

352 Furthermore, we analyzed livers for changes in pathology, and observed a reduction in the size 

353 of the lesions, hepatocyte necrosis, and inflammation in AM2 infected mice compared to wild- 

354 type infected (Fig. 5D). These data indicate that AM2 is attenuated compared to wild-type MF1V. 

355 To determine if primary infection with AM2 mutant virus protected mice from challenge 

356 with wild-type virus, we infected AM2 immunized and age-matched naive C57BL/6 mice 

357 intracranially with 6000 PFU of wild-type virus at 9 weeks post primary infection. We monitored 

358 body weight loss over the course of infection and found that naive mice lost significantly more 

359 weight starting at day 2 post-infection (Fig. 5E). In contrast, AM2 immunized mice did not lose 

360 weight and did not exhibit other symptoms of disease upon challenge. 13-week old C57BL/6 are 
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361 much more resistant to infection than 4-week mice, so as expected, no wild-type or AM2 

362 infected mice succumbed to the infection. These results indicate that the immunization with 

363 AM2 generated a sufficient immune response to protect the mice from subsequent wild-type 

364 challenge. 

365 

366 Discussion 

367 Modulating viral enzyme activity is a powerful approach for generating attenuated 

368 viruses. Coronaviruses encode a very large replicase polyprotein that contains multiple 

369 enzymatic activities such as proteases, ADP-ribose-1 ’’-phosphatase (ADRP), methyl-transferase, 

370 exonuclease, and RNA-dependent RNA polymerase (see Figure 1). Previous studies have shown 

371 that inactivation of the coronavirus papain-like protease by mutagenesis of the required catalytic 

372 cysteine residue or inhibition with a protease inhibitor, blocks viral replication (24, 37, 43). 

373 These approaches demonstrate the requirement for protease activity, but represent “all or 

374 nothing” approaches and do not allow for the study of reduced or impaired protease activity. 

375 Here, we describe an approach involving mutation of the domain adjacent to the papain-like 

376 protease domain, termed the ubiquitin-like domain (Ubl). Structural studies revealed that this 

377 domain is highly conserved in coronaviruses, although the function of this domain is unknown. 

378 We found that mutation of the Ubl domain destabilizes the adjacent protease domain, leading to 

379 loss of enzymatic activity over time (half-life ~30 min at 30°C). When this mutation is 

380 introduced into the virus, the mutant virus replicates efficiently in cell culture at 37°C, but is 

381 temperature sensitive and unable to replicate at high temperature (40.5°C). Indeed, we found 

382 that this virus is highly attenuated in infected mice but does replicate sufficiently to induce a 
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383 protective immune response against wild-type virus. Overall, our studies revealed a new domain 

384 that can be exploited for attenuating the pathogenesis of a coronavirus. 

385 Ubl domains have been identified as adjacent to many viral and cellular proteases, 

386 particularly the ubiquitin specific proteases (USPs), a family of deubiquitinating enzymes 

387 (DUBs) (44). Ubl domains can be present at the C-terminal or N-terminal end of USP catalytic 

388 domains, or even can be inserted into the catalytic domain. In addition, USPs may contain one 

389 or multiple Ubl domains (45). The Ubl domain of USP 14 is required for targeting this cellular 

390 DUB to the proteasome. At the proteasome, USP 14 removes ubiquitin molecules from proteins 

391 that are modified with K48-linked ubiquitin, which serves as a signal for proteosomal 

392 degradation (46). Alternatively, the Ubl may regulate function of DUB activity by interacting 

393 with other cellular proteins such as GMP synthetase, which enhances the activity of USP7 (47). 

394 Previous studies have shown that the presence of the Ubl domain of SARS-CoV PLpro is 

395 necessary for efficient expression and purification of the enzyme (37). To understand the 

396 mechanism by which mutations in the Ubl domain of MHV PLP2 affected function, we mapped 

397 the mutated residues V785 and V787 into the recently solved crystal structure of MHV PLP2 

398 (Fig. 6) (Chen et al., in preparation). Interestingly, the overall fold of the protease and Ubl-2 

399 domain is conserved between SARS-CoV PLpro and MHV PLP2. Valine 787, which is the 

400 focus of this study, is conserved between MHV and highly pathogenic coronaviruses such as 

401 SARS-CoV and MERS-CoV (Fig. IB). Valine 787 resides within a P-sheet that helps to 

402 maintain the structure of the Ubl. The side chains of valine 785 and valine 787 participate in the 

403 formation of a hydrophobic core of the Ubl, which stabilizes this domain. By mutating valine 

404 787 to serine, we introduced a polar side chain that is energetically unable to pack into the 

405 hydrophobic core because of significant desolvation penalties. Therefore, mutation of valine- 
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406 787 or valine 785 likely disrupts the structure of the Ubl domain thereby leading to the overall 

407 destabilization of PLP2 catalytic domain at elevated temperature. As a consequence, the Ubl 

408 mutant enzyme has decreased thermal stability that results in a loss of the specific activity of the 

409 enzyme over time in vitro. Furthermore, the destabilization of Ubl-2-PLP2 has an important 

410 impact on virus pathogenesis since we found that MHV-AM2 is attenuated in vivo. 

411 Interestingly, our studies revealed that the most dramatic reduction in AM2 viral titer in 

412 the mice was in the liver, not the brain, which was the site of inoculation. Since MHV-A59 is 

413 hepatotropic, it rapidly disseminates to the liver of infected animals (19). The reduction in viral 

414 titer in the liver and the decreased size of lesions in the liver (Fig. 5) could be due to either 

415 defective viral replication or a more rapid immune response to clear the virus from the liver. 

416 Similar phenotype of virus attenuated in the liver but not the brain was described for MHV-A59 

417 ns2 mutant (48, 49). In addition, a recent study by Zhang and colleagues showed that an nspl 

418 mutant virus is attenuated in the liver but replicates efficiently in the brain (50). The attenuated 

419 phenotype of the nspl mutant virus is hypothesized to be due to a more robust immune response 

420 to the virus. 

421 Recent studies on MERS-CoV showed that PLpro has DUB activity and that this activity 

422 is important for antagonism of innate immune response in transfected cells (31, 51). In addition, 

423 the arterivirus equine arteritis virus (EAV) PLP DUB activity has been shown to be involved in 

424 inhibition of innate immune response in infected cells suggesting that PLP DUB activity might 

425 be important for vims pathogenesis in vivo (52). Further studies are needed to determine the role 

426 of DUB activity in viral replication and pathogenesis and to determine if the attenuation of AM2 

427 is primarily immune mediated or due to defective replication and dissemination of the vims. 
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428 Coronaviruses encode multiple novel enzymatic activities and previous studies support 

429 the concept of modifying/mutating nonstructural proteins to generate attenuated viruses and to 

430 identify the role of the protein in virus replication and in causing disease in mice. Mutation of the 

431 N-terminal replicase product nspl in MHV results in an attenuated infection in mice (53, 54). In 

432 addition, mutation of nspl in SARS-CoV results in a decrease in antiviral signaling and virus 

433 titer in infected cells (55). Mutation of the catalytic residues of the ADP-ribose-1”-phosphatase 

434 (ADRP) domain does not diminish virus replication in mice but reduces production of the 

435 cytokine IL-6, an important pro-inflammatory molecule (56-58). Reduction in inflammatory 

436 response upon infection with ADRP mutant may contribute to decreased liver pathology (56). 

437 Genetically inactivating the enzyme ExoN, which has proofreading function, results in a 

438 hypermutation phenotype virus that is highly attenuated, even in aged, immunocompromised 

439 animals (59). Mutation of the 2’-0-methyl-transferase activity in nspl6 revealed that this 

440 activity is required to evade viral mRNA detection by the pattern recognition receptor MDA5 

441 (60, 61). Further, deletion or introduction of an inactivating mutation in nonstructural protein 2 

442 (ns2) in MFIV results in the activation of the RNase L pathway, which induces of an antiviral 

443 state in liver macrophages (62). All of the aforementioned studies reveal the critical roles of viral 

444 nonstructural proteins in virus replication and escape from immune surveillance. Our study 

445 reveals a role for the Ubl-2 domain in stabilization of the PLP2 catalytic domain of MFIV PLP2 

446 and provides a new target for viral attenuation and potentially for antiviral drug development. 

447 
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656 Figure legends 

657 Figure 1. MFIV-A59 Ubl-2 domain influences papain-like protease activity. A) Schematic 

658 diagram of MHV-A59 ORFs and the Ubl-2 and papain-like protease (PLP) domains within 

659 nonstmctural protein 3. Expression plasmids PLP2, and the VDVL region (aa 785-788) are 

660 indicated. B) Alignment of the predicted ubiquitin-like domains from MHV and MERS-CoV 

661 with the structural information from the Ubl domain of SARS-CoV (PDB 2FE8). C) F1EK293T 

662 cells were transfected with plasmids expressing PLP2 wild-type, PLP2 catalytic mutant CA, 

663 VDV Ubl-2 domain mutants in the presence of plasmid expressing the nsp2/3-GFP substrate. At 

664 24 hours post-transfection cells were lysed and analyzed by Western blot. D) PLP2 activity in 

665 live-cell assay. FIEK293T cells were transfected with pGlo-RLKGG and indicated PLP2 

666 expression plasmids. At 14 hours post infection GloSensor was added and luminesce was 

667 assessed hourly. E) FIEK293T cells were transfected with Flag-Ub expression plasmid, and wild- 

668 type (PLP2-WT) or indicated Ubl-2 mutants. Cells were lysed 24 hours post-transfection and 

669 analyzed by Western blot. The figure shows representative data from at least two independent 

670 experiments. 

671 Figure 2. Murine coronavirus AM2 is temperature sensitive for replication. A) Summary of 

672 the characteristics of Ubl mutant viruses recovered by reverse genetics. B) Representative 

673 plaques generated by wild-type MF1V (WT) or AM2 during incubation at the indicated 

674 temperature. DBT cells were inoculated with either WT or AM2, incubated at the indicated 

675 temperature and fixed and stained 48 hours post infection. Diameter of each plaque was 

676 measured and average diameter and standard deviation indicated for each virus; *p<0.05, 

677 **p<0.001. Figure shows representative data from at least two independent experiments. 
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678 Figure 3. Replication kinetics of AM2. . A) AM2 replication is similar to MHV-WT at 37°C. 

679 (A) DBT cells were infected with AM2 or WT at MOI 0.1 or (B) Bone marrow derived 

680 macrophages were infected with AM2 or WT at MOI 1. At indicated time points, supernatant 

681 was harvested and vims titer determined by plaque assay at 37°C. C) Temperature shift reduces 

682 yield of AM2. DBT cells were infected at the MOI 0.1 at 37°C with AM2 or WT. At the 

683 indicated time point post infection, cells were moved to 39.5°C until 14 hrs pi when supernatant 

684 was harvested and virus titer was determined by plaque assay at 37°C. The error bars represent 

685 standard deviation within single experiment. Figure shows representative data from at least two 

686 independent experiments. 

687 Figure 4. Ubl-2 domain is important for MFIV PLP2 thermal stability. A) SDS-PAGE 

688 analysis of the final purified PLP2-V787S enzyme. Wild-type and PLP2-V787S mutant PLP2 

689 were expressed and purified from E. coli as described in Materials and Methods. B) 

690 Temperature-dependent inactivation of wild-type and PLP2-V787S mutant. The specific activity 

691 of wild-type (WT) and the V787S mutant PLP2 enzymes were measured after incubation at 25°C 

692 and 30°C for different time periods. The specific activities were then normalized to the activity 

693 at 0 min (Rate/Rateo: rate at time t over initial rate). Kinetic data for the V787S mutant 

694 incubated at 30°C were fit to a first-order exponential decay model (Rate t /Rateo = e" kt ). The 

695 calculated inactivation rate constant ki„ act for the V787S mutant at 30°C is 0.025 ± 0.001 min" 1 

696 and the half-life (ti/?) is 27.7 ± 1.1 min. C) CD melting curves of WT and PLP2-V787S mutant 

697 PLP2. The thermal stability of WT and V787S mutant PLP2 was determined by measuring the 

698 CD signal at 220 nm as a function of temperature at a step interval of 0.4°C and at a rate 

699 0.5°C/min. Three independent experiments were performed for both WT PLP2 (grey) and 

700 V787S mutant (black). 
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701 Figure 5. MHV AM2 is attenuated and generated protective immunity in mice. A) C57BL/6 

702 mice were infected with 600 pfu WT icMHV or AM2 mutant intracranially and monitored for 

703 survival over time (N=7 for each group). B) Virus titers in brains of mice infected with 600 pin 

704 WT icMHV or AM2 were determined at indicated time points p.i. by plaque assay (N=4 to 5). 

705 Error bars represent standard deviation. C) RNA levels for viral replicase gene and cellular IL-6 

706 were determined using qRT-PCR in livers at day 5 post infection (N=6). RNA levels were 

707 normalized to actin and amount of transcript in WT infected mice was set to 1. Error bars 

708 represent SEM. Statistical analysis performed using Student’s t test. D) Representative samples 

709 of hematoxylin and eosin staining of formalin-fixed liver samples at day 5 post infection from 

710 mice infected with WT icMHV or AM2 (N=6 for each group, magnification X10). E) C57BL/6 

711 mice immunized with AM2 mutant and naive age-mached controls were challenged with 6000 

712 pfu icMHV intracranially 9 weeks post primary infection. The mice were monitored for body 

713 weight loss. Error bars represent SEM.****, p<0.0001, two-way ANOVA. 

714 Figure 6. X-ray structures of the MHV and SARS-CoV Ubl-2 domains. A superposition of 

715 the Ubl-2 domains from MHV PLP2 (blue) and SARS-CoV PLpro (orange, PDB 2FE8) is 

716 shown as a stereoview and resulted in RMSD of 0.74 A. Side chains of MHV Ubl-2 residues 

717 V785 and V787, and their corresponding residues in SARS-CoV are shown as sticks. 


D 

o 


<D 

C 



A) MHV 



B) 


SARS-CoV 

SARS-CoV 

MHV 

MERS-CoV 


TTQQLTll 


P3 &4 



111 


P5 


|\/T KIXPHVNH K0K TFF7LPS0 
KVRCCXI YK0K7FFQYCDL 

ot i pdekqnBhIslIyIladnl 


C) 


D) 


o 

o 

=3 

■O 


2 

o 


20 

16 

12 


4 

0 



E) 



Hours GloSensor Incubation 



PLP2 

Calnexin 


Ubl-2 



Figure 1. MHV-A59 Ubl-2 domain influences papain-like protease activity. A) Schematic diagram of MHV-A59 
ORFs and the Ubl-2 and papain-like protease (PLP) domains within nonstructural protein 3. Expression plasmids 
PLP2, and the VDVL region (aa 785-788) are indicated. B) Alignment of the predicted ubiquitin-like domains from 
MHV and MERS-CoV with the structural information from the Ubl domain of SARS-CoV (PDB 2FE8). C) HEK293T 
cells were transfected with plasmids expressing PLP2 wild-type, PLP2 catalytic mutant CA, VDV Ubl-2 domain 
mutants in the presence of plasmid expressing the nsp2/3-GFP substrate. At 24 hours post-transfection cells were 
lysed and analyzed by Western blot. D) PLP2 activity in live-cell assay. HEK293T cells were transfected with pGlo- 
RLKGG and indicated PLP2 expression plasmids. At 14 hours post infection GloSensor was added and luminesce 
was assessed hourly. E) HEK293T cells were transfected with Flag-Ub expression plasmid, and wild-type (PLP2- 
WT) or indicated Ubl-2 mutants. Cells were lysed 24 hours post-transfection and analyzed by Western blot. The 
figure shows representative data from at least two independent experiments. 
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Virus 

Ubl Mutation 

Titer of recovered virus 
(pfu/ml) 

Plaque size 

WT 

none 

10 7 

regular 

AMI 

VDVL/SSSS 

10 3 

pin point 

AM2 

V787S 

10 7 

small 

AM3 

V785S 

10 7 

small 


B) 

WT 

AM2 


37.0°C 



3.98mm ± 0.58 — 

M 

3.26mm + 0.40 — 


39.5°C 40.5°C 



2.46mm+ 0.11 0.00mm ± 0.00 


Figure 2. Murine coronavirus AM2 is temperature sensitive for replication. A) Summary of 
the characteristics of Ubl mutant viruses recovered by reverse genetics. B) Representative 
plaques generated by wild-type MHV (WT) or AM2 during incubation at the indicated 
temperature. DBT cells were inoculated with either WT or AM2, incubated at the indicated 
temperature and fixed and stained 48 hours post infection. Diameter of each plaque was 
measured and average diameter and standard deviation indicated for each virus; *p<0.05, 
**p<0.001. Figure shows representative data from at least two independent experiments. 
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Figure 3. Replication kinetics of AM2. A) AM2 replication is similar 
to MHV-WT at 37°C. (A) DBT cells were infected with AM2 or WT at 
MOI 0.1 or (B) Bone marrow derived macrophages were infected 
with AM2 or WT at MOI 1. At indicated time points, supernatant was 
harvested and virus titer determined by plaque assay at 37°C. C) 
Temperature shift reduces yield of AM2. DBT cells were infected at 
the MOI 0.1 at 37°C with AM2 or WT. At the indicated time point 
post infection, cells were moved to 39.5°C until 14 hrs pi when 
supernatant was harvested and virus titer was determined by plaque 
assay at 37°C. The error bars represent standard deviation within 
single experiment. Figure shows representative data from at least 
two independent experiments. 
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Figure 4. Ubl-2 domain is important for MHV PLP2 thermal stability A) SDS-PAGE analysis of 
the final purified PLP2-V787S enzyme. Wild-type and PLP2-V787S mutant PLP2 were expressed 
and purified from E. coli as described in Materials and Methods. B) Temperature-dependent 
inactivation of wild-type and PLP2-V787S mutant. The specific activity of wild-type (WT) and the 
V787S mutant PLP2 enzymes were measured after incubation at 25°C and 30°C for different time 
periods. The specific activities were then normalized to the activity at 0 min (Rate,/Rate 0 : rate at 
time t over initial rate). Kinetic data for the V787S mutant incubated at 30°C were fit to a first-order 
exponential decay model (Rate,/Rate 0 = e _kt ). The calculated inactivation rate constant k inac , for the 
V787S mutant at 30°C is 0.025 ± 0.001 min' 1 and the half-life (t 1/2 ) is 27.7 ± 1.1 min. C) CD melting 
curves of WT and PLP2-V787S mutant PLP2. The thermal stability of WT and V787S mutant PLP2 
was determined by measuring the CD signal at 220 nm as a function of temperature at a step 
interval of 0.4°C and at a rate 0.5°C/min. Three independent experiments were performed for both 
WT PLP2 (grey) and V787S mutant (black). 
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Figure 5. MHV AM2 is attenuated and generated protective immunity in mice A) C57BL/6 mice 
were infected with 600 pfu WT icMHV or AM2 mutant intracranially and monitored for survival over 
time (N=7 for each group). B) Virus titers in brains of mice infected with 600 pfu WT icMHV or AM2 
were determined at indicated time points p.i. by plaque assay (N=4 to 5). Error bars represent 
standard deviation. C) RNA levels for viral replicase gene and cellular IL-6 were determined using 
qRT-PCR in livers at day 5 post infection (N=6). RNA levels were normalized to actin and amount of 
transcript in WT infected mice was set to 1. Error bars represent SEM. Statistical analysis performed 
using Student’s f test. D) Representative samples of hematoxylin and eosin staining of formalin-fixed 
liver samples at day 5 post infection from mice infected with WT icMHV or AM2 (N=6 for each group, 
magnification XI0). E) C57BL/6 mice immunized with AM2 mutant and naive age-mached controls 
were challenged with 6000 pfu icMHV intracranially 9 weeks post primary infection. The mice were 
monitored for body weight loss. Error bars represent SEM.****, p<0.0001, two-way ANOVA. 
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Figure 6. X-ray structures of the MHV and SARS-CoV Ubl-2 domains. A superposition of the 
Ubl-2 domains from MHV PLP2 (blue) and SARS-CoV PLpro (orange, PDB 2FE8) is shown as a 
stereoview and resulted in RMSD of 0.74 A. Side chains of MHV Ubl-2 residues V785 and V787, 
and their corresponding residues in SARS-CoV are shown as sticks. 
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Table 1. Primers 


Forward Primer 

Reverse Primer 

PLP2-CA 5' G TCA AAC AAC AAC GCC TAC ATC AAC G 

5' CGT TGA TGT AGG CGT TGT TGT TTG AC 

SSSS 5' GGC CAA TAA GAG CAG TAG TAG TTG CAC CGT CGA CGG 

5' CCG TCG ACG GTG CAA CTA CTA CTG CTC TTA TTG GCC 

V785S 5' GGC CAA TAA G AGT GAT GTG CTG TG CAC CGT CGA CG 

5' CGT CGA CGG TGC ACA GCA CAT CAC TCT TAT TGG CC 

D786K 5’ GGC CAA TAA G GTC AAA GTG CTG TG CAC CGT CGA CG 

5' CGT CGA CGG TGC ACA GCA CTT TGA CCT TAT TGG CC 

V787S 5' GGC CAA TAA G GTC GAT TCG CTG TG CAC CGT CGA CG 

5' CGT CGA CGG TGC ACA GCG AAT CGA CCT TAT TGG CC 


AMI 5’ CTC GCT AAT AAG AGT TCG TCC TCG TGT ACT GTT GAT GG 5’ CCA TCA ACA GTA CAC GAG GAC GAA CTC TTA TTA GCG AG 

AM2 5’ CTC GCT AAT AAG GTT GAT TCC TTG TGT ACT GTT GAT GG 5’ CCA TCA ACA GTA CAC AAG GAA TCA ACC TTA TTA GCG AG 

AM3 5’ CTC GCT AAT AAG TCT GAT GTC TTG TGT ACT GTT GAT GG 5’ CCA TCA ACA GTA CAC AAG ACA TCA GAC TTA TTA GCG AG 

Replicase 5’ GAA TGC CCG GGC GGG ATT TTT GTA TCC 5’ CCA CAA GAT CTG CCT CGG ACA AAT C 






